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ABSTRACT 



Context. Hard X-Ray bremsstrahlung continuum spectra, such as from solar flares, are commonly described in terms of 
power-law fits, either to the photon spectra themselves or to the electron spectra responsible for them. In applications 
various approximate relations between electron and photon spectral indices are often used for energies both above and 
below electron low-energy cutoffs. 

Aims. We examine the form of the exact relationships in various situations, and for various cross-sections, showing that 
empirical relations sometimes used can be highly misleading and consider how to improve fitting procedures. 
Methods. We obtain expressions for photon spectra from single, double and truncated power-law electron spectra for a 
variety of cross-sections and for the thin and thick target models and simple analytic expressions for the Bethe-Heitler 
cases. 

Results. We show that above a low-energy cutoff the Kramers and Bethe-Heitler results match reasonably well with 
results for exact cross-sections up to energies around 100 keV; that below the low-energy cutoff, Kramers and other 
constant spectral index forms commonly used are very poor approximations to accurate results; but that our analytical 
forms are a very good match. 

Conclusions. Analytical forms of the Bethe-Heitler photon spectra from general power-law electron spectra are an 
excellent match to exact results for both thin and thick targets and they enable much faster spectral fitting than 
evaluation of the full spectral integrations. 

Key words. X-Ray; spectrum; power-law; spectral index; bremsstrahlung; cross-section; Sun: corona-Sun; activity-Sun: 
flares-Sun: X-rays, spectra 



1. Introduction 

Hard X-ray bremsstrahlung spectra are important diagnos- 
tics of flare electron acceleration and propagation - e.g., 



Brownl(ll971h: Lin & Schwartz! (11987 
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enabled b y the high resolution spectra being observed by 
RHESSI (|Lin et al.l 12002.) . which handles very large spec- 



tral and dynamic ranges. At low energies (a few keV for mi- 
croevents and up to 20 or so keV for large flares), the spec- 
trum i s usually consistent w ith isothermal bresmsstrahlung 
- e.g. (|Holman et al.l [20031 ) - while at higher energies it is 
usually consistent with bremsstrahlung from (a sometimes 
broken) power-law electron spectrum with a low-energy 
cutoff. One should bear in mind that other, broadly sim- 
ilar, forms (e.g. shifted power-laws {E + E„)~"') are also 
consistent with the data and that the data /(e) to be con- 
sidered are, in general, those after appli cation of corrections 
for the albedo spectrum contribution (jKontar et al.l 120061 ) 
as well as instrumental effects. However, here we focus on 
the properties of truncated power-law fits (including broken 
power-laws) since these are so widely used. 
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We consider here mainly the properties of the higher en- 
ergy component and discuss the relationship between the 
'local' spectral indices 6{E) of the source electrons and 7(e) 
of th e observed photons. In data analysis it is quite common 
fe.g. iHannah et al]|2007l ) to assume constant values of 7(e) 
and of S(E) over specific finite energy ranges and definite 
linear relationships between 7 and S in these - for exam- 
ple, in the energy ranges below and above the electron cut- 
off. However, as we show below, for general bremsstrahlung 
cross-sections, most such relationships are at best approxi- 
mate and the various ad hoc relationships used can be quite 
misleading. Our aim is to show the exact form of these rela- 
tionships for various cross-sections and derive analytic ex- 
pressions simple enough for easy use in fast spectral analysis 
software. 

For both the thin- and thick-target models this is ac- 
complished by first showing numerically that the Bethe- 
Heitler approximation for the bremsstrahlung cross-section 
is reliable enough to reproduce accurately the true pho- 
ton spectrum corresponding to a truncated power-law elec- 
tron spectrum. Then for both models we obtain exact an- 
alytical expressions for /(e) and 7(e) based on the Bethe- 
Heitler cross-section. Finally these expressions are used to 
best fit simulated and measured photon spectra and de- 
termine electron spectrum parameters. The effectiveness of 
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this method is assessed by comparisons with two differ- 
ent fitting approaches currently employed: (1) fitting with 
a numerical ex pression for /(e) ob tained by full numeri- 
cal integration (jHolman et al.ll2003D of truncated power-law 
electron spectra F{E) and J^o{Eo) (thin and thick target, 
respectively); (2) fitting data on /(e) with a broken power- 
law characterized by two distinct constant photon spec- 
tral indices, one abo ve and one below a 'knee' energy (e.g. 
I Hannah et all [20071 ) . With respect to the first approach, 
we find that our new method yields equally good values 
but with substantially higher computational speed (much 
higher in the case of thick targets). With respect to the 
second approach, our method is no faster but much more 
accurate and meaningful. In fact the /(e) used in approach 
(2) is unphysical, corresponding to no real F{E) so it may 
give excessive values due to large residuals near the knee. 

Section 2 provides the general equations for 
bremsstrahlung spectra for both coUisionally thin and 
thick target sources. Section 3 establishes our notation 
for a power-law F{E) of constant S truncated below the 
low-energy cutoff Ei^ and shows how arbitrarily trun- 
cated and multiple (broken) power-laws can be expressed 
in terms of these. Section 4 defines the Kramers and 
Bethe-Heitler approximate bremsstrahlung cross-sections 
and contrasts the results they give for /(e) and 7(e) for 
power-law electron spectra compared with that for the 
exact cross-section. In Sections 5 and 6 we obtain analytic 
expressions, in both thin- and thick-target cases, for /(e) 
and 7(e). In Section 7 we report numerical tests of the 
speed and accuracy of using these to fit real and simulated 
data, as compared with other approximate methods and 
with full integration. Section 8 summarises our conclusions. 

2. Thin and thick target bremsstrahlung and 
energy losses 

For a general inhomogeneous optically thin source 
of plasma density n(r) and electron flux energy spectrum 
F{E, r) in volume V, the bremsstrahlung photon flux en- 
ergy spectrum /(e) (cm~^s~^ p er unit e at Earth distance 
R) can be written (jBrownlll97lh : 



He) 



with 



nV 
47r/?2 



F{E)Q{e,E)dE 



n= / ndV/V 
'v 

and 



F{E) 



nir)F{E,r)dV/{nV) 



(1) 



(2) 



(3) 



where Q{e,E) is the bremsstrahlung_ cross-section differen- 
tial in photon energy e. In general, F{E) and Q(e, E) have 
to be treated as anisotro pic and Equation (HI) involves an 
integral over solid angle (|Brownlll97^lMassone et al.ll2004D 
though most data treatments assume source isotropy. 

In a purely coUisional thick target, F{E) is related to 
the thick-target injection rate spectrum J-o{Eq) (electrons 
per second per unit injection energy Eq) through equation 



KnV Je 



(4) 



regardless of Q(e, E), with K = 27r e^A and A the Coulomb 
logarithm - iBrown fc Emsli^ (|1988f ). Though they consider 
only the case K = l-ne^k for collisional losses only in a uni- 
formly ionized target. Equation applies to any energy 
loss rate coefficient K{E) if the loss rate can be written in 
the form dE/dN = —K{E)/E where N is the column den- 
sity along the electron path such as for collisional losses at 
high E with relativistic correction. Synchrotron losses can- 
not be written in this way unless the magnetic field and elec- 
tron pitch angle distribution do not vary along the path, but 
these also only matter at high energies. The most serious 
approximation involved in using constant K for lower ener- 
gies is in (common) neglect of the fact that K varies with 
target hydrogen ionization x ([Brownl 119731 : lEmslid Il978f ) 
being K = 27re''A(l -I- ax) /{1 + a) with a ~ 1.6. This is im- 
portant around the energies of electrons /?, w {2KNf,)^/^ 
stopping around the flare transition zone depth N^,. Apart 
from early stages of flares, prior to much evaporation, /?, 
is well above typical values of Ei considered here. Since we 
are concerned mainly with small E around low-energy cut- 
off values, henceforth we address only the case of constant 
K and take A = 25. 



3. Single and broken power-laws 

We consider first the widely used single power-law form 
with low-energy cutoff for the thin target F{E) 



F{E) 




(5) 



where Fi = F{E)dE is the total mean electron flux at 

E > El, a l ow-energy cutofl0i and d is the (thin target; 
iBrownl ( 1971[ )) constant electron spectral index. 

Before considering the photon spectral properties of 
such single power-law electron spectra we note that our 
results for these can easily be generalised to fitting of dou- 
ble (broken) power-laws in F{E) with lower and upper cut- 
off energies. The following decomposition expressions apply 
equally well to any broken and truncated power-laws such 
as in photon space /(e). The general case is 



FiE) 





AE-^^ 
AE, 




'E 



-S2 



E<E, 
Ec<E <Eb 

Eb<E<Ea 
E>Ea 



(6) 



where A, Ec, Eb, Ea, Si, 62 are constants. Writing the 
parametrized single power-law Fpi{c, d, E*) as 



Fpiic d, E*] 



cE--^ E>E* 
otherwise , 



(7) 



(8) 



F{E) in ^ can be always written as 

F{E) = Fpi{A,Si,E,)-Fpi{A,6i,Eb) + 

Fpi{AE-''+'\62, Eb) ~ Fpi{AE-''+'\62, Ea). 

Hence the corresponding /(e) can be found simply as the 
sum and difference of the relevant /(e) expressions for single 
power-laws. 



^ If one wishes to use an electron flux reference energy E, 
distinct from the cutoff energy Ei one must replace Fi in Eq. 
® by F.(S. /El 



Brown et. al.: Fast Spectral Fitting 



3 



For the thick target model, we have to revisit the prob- 
lem and evaluate the form of F{E) and hence of /(e) for 
a truncated power law in To{Eq), not in F{E). Note also 
that we have to distinguish between the spectral index 5o 
for a pure power-law injection spectrum J^q{Eq) from the 
index 5 for F{E). For an injection spectrum truncated at 
Eq < Eoi 




Eo > Eqi 
Eq < Eqi 



(9) 



where J^oi — J ^{Eo)dEQ and equation (jH) gives 



F{E)^ 



KnV 



Eoi 
E 
Eoi 



-Sa+2 



E > Eoi 
E < Eqi 



(10) 



i.e., the relation between and S at E > Eqi being d — 
6q — 2. To find F{E) from a broken and truncated power- 
law form of the injected thick target To{Eq), one would first 
use the analogy of expression ([8]) to get the total J^o{Eo) 
as the sum of a set of single power- law J^q{Eq) forms and 
use expression j[TU)) for each term in the sum to get the 
corresponding F{E). 



4. Cross-sections 

The bremsstrahlung cross-section Q{e, E) can be writ- 
ten as 



eE 



(11) 



incorporating a high Z element correction factor AzZ^ 
in Qo, the Gaunt factor g(e, E) depending on the actual 
cross-section used. For the Kramers approximation 

q{e,E)^qK{e,E)^l- (12) 

while for the non-relativistic Bethe-Heitler approximation 



q{e,E) = qBHie,E) = log 



1 + y/l-e/g 
1 - VI - ^/E 



(13) 



The most widely used isotropic formula (neglect- 
ing electron-electron bremsstrahlur i g wh ich is impor- 
tant at high energies (jKontar et al.l l2007h ) is q{e,E) — 
Q:^ rn(^-,E) correspo nding to the 3BN formula from 
iKoch fc Mot^ (I1959D used in inversion and forward fits 
in iPiana et"al] (l2003h: iKontar et all (|2005f ): iBrown et all 
(|2006f) : iJohns k LinI (|l992f) . In the case of a truncated 
power-law, the sensitivity of the predicted thin- and thick- 
target photon spectra to the form of the cross-section used 
is illustrated in Figure 1 by numerically computing 7(e) for 
6 = 6 — b] 5q — b, 5q = 7 and various q. It is clear firstly 
that, particularly below the cutoff, Bethe-Heitler results are 
a rather good approximation for computing local photon 
spectral indices 7(e). Second, above and especially below 
the cutoff 7(e) is not in any way constant, as sometimes 
assumed ()Hannah et al. I [20071 ). Figure 1 thus clearly sug- 
gests that the analytical expressions of J(e) and 7(e) based 
on the Bethe-Heitler formula could be useful for spectral 
computation in forward spectral fitting. 

Note that Figure 1 shows results for only one value of Ei 
(Eqi) and presents them only as functions of e/Ei (e/Eoi) 



whereas one might in general expect results to depend on e 
and El (Eqi) separately. However, we carried out a range 
of test calculations for several different Ei {Eqi ) in the few 
deka-keV range and found that results were, to a very good 
approximation, functions only of e/Ei (e/i?oi). Secondly, 
the most general Q is actually anisotropic and the above 
expressions have to be ge neralised to integrat e over electron 
angle as well as energy ([Massone et al.|[2004[ ). However, the 
effect of this on J(e) is small at low energies and in any 
case is mainly an e-independent scaling rather than a spec- 
tral effect. In fact some numerical experiments showed that 
there is little effect on our conclusions of using anisotropic 
F,Q. 



5. Thin target spectra and spectral indices 

Using Equation pT|) . Equation ^ becomes 

Qquic^ nV 



FiE)q{e,E) 



dE 

"e 



(14) 



We are interested here in results for photon spectra /(e) 
when F{E) is a single power-law with constant 5 and a low- 
energy cutoff (Equation ([5])), for various forms of q{e.,E) 
considered, and also on the form of local photon spectral 
index defined by (Brown and Emslie 1988) 



7(e) 



e dl d log / 



I de 



dloge 



(15) 



for comparison with the use (e.g. iHannah et al.l l2007f ) of 
constant 7 approximations. For any q{e,E) equation ([14]) 
can be written 



e 



minllXEi/e)"] 



j{e,Ei/x^/^)dx (16) 



where x — [Ei/EY , and 



C 



Qqitic^ nVFi 
47r/?2 eI~ 



(17) 



Note that in special cases where q{e, E) — q{e/E) only, 
7(e) takes the form 



7(e) 




log Ulq{x^/'elEi) 



e>Ei 
dx e < El 



(18) 



5.1. Kramers cross-section 

For Kramers q{e, /?) — 1 and we have immediately 

/(e) = i-ic^(m'^>^i 

^ ' S ell e<Ei 



and 

7(e) 



S+1 e>Ei 
1 e<Ei . 



(19) 



(20) 
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0.1 1.0 10.0 0.1 1.0 10.0 

(c) (d) 



Fig. 1. Local bremsstrahlung photon spectral index 7(e) behavior computed numerically for three different q{e,E) - 
Kramers, Bethe-Heitler and 3BN. (a) Thin target with (5 = 3. (b) Thin target with 5 = 5. (c) Thick target with Sq = 5. 
(d) Thick target with So = 7. 



5.2. Bethe-Heitler cross-section 

The spectrum and the spectral index can be written 
solely in terms of the electron spectral index S and the 
dimensionless parameter a = e/Ei. Integration by parts 
leads to: 



I{S,a) = 



-C 



a > 1 

(21) 
a < 1 



with B{a,f3) = - xf-^dx and B^{a,(i) 

t"~^{l — t)^~^dt the incomplete beta function. 
For the spectral index: 



6. Thick target spectra and spectral indices 

Inserting into (HH) leads to 

/J ^ ^ q{e,Eoi/y^o~-)dy e>Eoi 

J^q{e,Eoi/y^)dy+ 
{do-2)l)^^^qie,iEoi)d^ e<Eoi 

where D = ^47^3$^; the integrals on the right hand side 
of ([23|) being expressed below in terms of the dimension- 
less parameter a = e/i?oi, the constant D and the spectral 
index Sq for the truncated electron power- law !Fq{Eq). 



log 



l + ^/T 



a> 1 6.1. Kramers cross-section 

a < 1 (22) 



Integration by parts with Kramers unity Gaunt factor 



We note that the previous analytical formulas for I{S, a) 
and 7(5, a) are written in terms of beta functions and in- 
complete beta functions. The computation of these func- 
tions is included in standard library routines for data visu- 
alization and analysis, making the implementation of the 
exact formulas (|2ip and easy and fast. 



gives 



1(60, a) 



and 



D ( a-''«+i a > 1 



l + (^0-2)(l~a) 



a < 1 



, So-1 a > 1 



(24) 



(25) 
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6.2. Bethe-Heitler cross-section 

Integration by parts with the Bcthe-Heitler Gaunt fac- 
tor qBH leads to 



r B(5o-2,i)[i] 



/(<5o,a) = 
and 

7(^0, a) = 
where 



D 1 



So-2 



6o-2-a\ [^o-l-^«]log 



So- I 

Ai+A2BJSo~2,^) 
A3-\/T^+A4B„((5o-2,i) 



a > 1 
a < 1, 



Al = Ai(5o,a) 
A2=A2{SQ,a) 

and 

A4 = A4{do,a) 



Sq-1 l + ^/l-a 
<5o-2 °^1-VT^ ' 

^0-1 1 
<5o-2a'5«-2 ' 

<5o-l 



So-2 2) °^1-VT~~^ 



1 



So -2 a' 



2 



a > 1 
+(26) 
a < I 

(27) 

(28) 
(29) 
(30) 

(31) 



Thus, as in the thin target case, integral expressions for 
I{So, a) and ^{Sq, a) can be written analytically in terms of 
beta functions and incomplete beta functions. 

7. Tests against data and comparison with other 
fitting methods 

To show the usefulness of our formulation in terms of accu- 
racy and speed, we have tried it out on simulated and real 
data. Our specific goals here for simulated noisy data are 
to determine : 

1. How good and how fast is use of our functional approxi- 
mations to /(e) in fitting data to estimate [S,nVFi, Ei] 
or [Jo, -^01 J -£'01 ] for a_[thin] or [thick] target single 
truncated power-law [F{E)], or [J^o{Eq)], in compari- 
son with other fitting routines. To answer this we use 
simulated data /(e) generated with NASA SolarSoft 
(SSW) routines using the full bremsstrahlung cross- 
section and numerical integ ration over E or Eo, respec- 
tively (|Holman et al.l 120031 ). 

2. How well do the speed and accuracy compare with 
other approaches? In particular with : (a) iHolman et al.l 
(|2003l ) who, in carrying out the best fit parameter 
searches, perform full integrations with the full cross- 
section in each iterative step; and (b) an approach 
which, instead of fitting the /(e) predicted for an elec- 
tron power law with cutoff, fit a parametric piecewise 
power-law photon /(e) with distinct constant photon 
spectral indices 71,72 at e <,> tb, wi th 71 either a 
free fit parameter or prescribed, e.g. as in I Hannah et al.l 
(|2007t ) where 71 — 1.5 (to be compared with 7(e) in 



Figure 1). Loosely speaking, for example in the case of 
thin targets, the value of the photon break energy e;, 
is meant to refiect an electron low-energy cutoff energy 
El, though in reality such a photon spectrum does not 
correspond to any real (non-negative) electron spectrum 
except for Kramers cross-section in which case 71 = 1, 
72 = J + 1- 

3. How each of the above findings changes when we add a 
reasonable isothermal contribution to /(e) for thermal 
parameters EM and T. 

We have carried out these simulated data comparisons 
for the following parameter sets: S,So — 3,5, Ei,Eoi = 
10,20,30 keV, nVFi = 5,20 x 10^^ electrons cm^^ g-i. 
Toi = 5, 20 X 10^5 electrons s'^, EM = 0.5, 1 x lO''^ cm'^, 
and kT = 1, 1.5 keV. 

/(e) were generated using SSW routines for thin and 
thick target (f_thin.pro and f_thick.pro, respectively) 
and optionally with isothermal component (f_vth.pro) in 
the 3 - 100 keV energy range with 1 keV energy bin- 
ning. Then, these /(e) were converted to counts using the 
RHESSI detector response matrix for attenuator state 0. 
Finally, we added Poisson noise. Such simulated count spec- 
tra were then fitted within the OSPEX enviroimient0 in the 
3-100 keV range. 

The results of comparisons for the simulated spectra 
without a thermal component are as follows: 

1. For both thin and thick cases, acceptable fits in terms of 
reduced x^j X^j and normalised residuals were found for 
all simulated spectra using expressions ([2T|) and ((26)l . re- 
spectively. The fitted parameters were close to the input 
ones - within ^ 10% or less. 

Concerning the computational time (see Table 1) for the 
fitting procedures, the SSW thin fit routine was about 2 
times slower than our method for similar accuracy while 
the SSW thick fit routine is 10-20 times slower than our 
method. 

2. As regards the matter of trying to get a meaningful 
fit to the actual form of /(e) from truncated power- 
laws F{E), !Fo{Eo) by using double power-law fits to 
the photon spectra (constant 7 above and below some 
break energy eb - bpow . pro) we found that this failed to 
produce an acceptable overall fit (generally, xt > 2) and 
that normalised residuals clustered near e ~ Ei,Eoi. 
Using a fixed value of 71 such as 1, 1.5, 1.7 did not help. 

Adding a plausible isothermal component with the 
above values for parameters EM and T to the thin and 
thick target spectra modifies the fit behaviour described 
above only if the thermal component contributes signifi- 
cantly to or dominates the spectrum at e ^ Ei,Eoi. For 
such spectra, e.g. thin-target case EM = 1 x lO''^ cm~^, 
kT = 1.5 keV, S ^ 3, Ei ^ 10 keV, the fitting functions 
introduced in this paper do give an acceptable fit but only 
an upper limit on Ei can be obtained. This limit is close to 
the energy where the thermal spectrum steepens and falls 
below the non-thermal part, fits with smaller Ei being also 
consistent with the data since lost in the dominant thermal 
emission. 

Photon spectra with a thermal component can also usu- 
ally be fitted with a nonthermal component close to a sin- 
gle power-law so the value of Cb obtained by double photon 
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SSW-thick Eq. (I26j 


SSW-thin 


Eq. ^ 


6, So = 3 








p TP —in 


7 4^ n 791 
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n 79fi 










p p — on 






n 79fi 


(5 <5n — 3 








TP fr _ on 


q QQ n 


u. yo -L 


n ^QQ 


^0 — 5 








£i,Soi = 10 


9.99 0.721 


1.36 


0.724 


S, So = 5 








Ei,Eoi = 20 


7.46 0.594 


1.07 


0.595 


(5, (5o = 5 








Bi.Soi = 30 


7.05 0.583 


1.20 


0.597 



Table 1. Computational cost test. The test was performed 
on a 3.2 GHz PC with 1GB RAM. Times are in seconds. 
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Fig. 2. Thin and thick target fits to February 20 2002 flare 
spectrum. The labels show the comparison between SSW 
and our approach. 



power law fits is not a good indicator of Ei, Eqi. Depending 
on the combination of thermal and non-thermal parts, ac- 
ceptable fits using bpow yield Cb which can be either lower or 
higher than the input Ei , Eqi . Therefo re, using eb for an es - 
timate of the non-thermal energy (e.g. lHannah et ani2007l ) 
can be misleading. On the other hand, our expressions gen- 
erally give El, Eqi much closer to the input Ei,Eqi than 
eb is. 

The proposed expressions have been also tested on two 
cases of real data. Figure 2 shows thin and thick fits to 
the 20-Feb-2002 11 UT flare and compares the parameters 
obtained by our and SSW expressions. Both thin and thick 
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Fig. 3. Low-energy cutoff and total non-thermal electron 
energy in 4s time intervals for the June 2 200 2 flare as 
derived from our approach, see also Figure 4 in ISui et all 
(I2007h . 

fits give similar F{E) and To{Eo). This flare was near the 
limb so required no albedo correction. 

Next, we applied our thick target expression to the early 
impulsive phase of the flare of 02-Jun-2002 which shows 
flattening and evidence of a low-energy cutoff at Ei above 
the t hermal component in the 18 - 38 keV range (jSui et al.l 
I2007D . In this case albedo correction was applied, as is es- 
sential for such events. Figure 3 shows the time evolution 
Ei{t) obtained for the time variation of the best fit low- 
energy cutoff and for the corresponding total non-thermal 
electron power (as total energy per 4 sec integrati on). These 
curves are closely comparable with those found bv lSui et al.l 
()200 7') in their Figure 4. 

Thin and thick target formulas (Equations (|21|26p ) 
have been incorporated into the SSW tree. Prospective 
users may access them as OSPEX fitting functions named 
photon_thin.pro and photon_thick.pro. 

8. Conclusions 

We have shown that results for thin and thick target 
bremsstrahlung photon spectra /(e) from power-law elec- 
tron spectra with constant index 6 and low-energy spec- 
tral cutoff El, Eqi obtained using the Bethe Heitler cross- 
section are very close to those from the exact cross-section, 
at photon energies both above and below the cutoff, at least 
at keV to deka-keV energies. We have shown further that 
the Bethe Heitler expressions allow the bremsstrahlung in- 
tegrals for /(e) to be written as analytic forms in terms of 
beta functions of S and a only, which are part of standard 
numerical packages, and that these give results very close 
to the exact /(e). 

For both the thin and especially for the thick target 
models, we find that this formulation enables equally accu- 
rate but much faster spectral fitting than evaluation of the 
full spectral integrations which will be valuable in analysis 
of bremsstrahlung HXR data such as from RHESSI. This 
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fast fit approach can also replace the use of fitting unphys- 
ical broken power-laws in /(e). 
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